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The metal content of planet hosting stars is an important ingredient which 
may affect the formation and evolution of planetary systems. Accurate stellar 
abundances require the determinations of reliable physical parameters, namely 
the effective temperature, surface gravity, microturbulent velocity, and metal- 
C/3 ■ licity. This work presents the homogeneous derivation of such parameters for 

a large sample of stars hosting planets (N=117), as well as a control sample of 
disk stars not known to harbor giant, closely orbiting planets (N=145). Stellar 
parameters and iron abundances are derived from an automated analysis tech- 
nique developed for this work. As previously found in the literature, the results 
in this study indicate that the metallicity distribution of planet hosting stars is 
more metal-rich by ~0.15 dex when compared to the control sample stars. A 
segregation of the sample according to planet mass indicates that the metallicity 
distribution of stars hosting only Neptunian-mass planets (with no Jovian-mass 
^ planets) tends to be more metal-poor in comparison with that obtained for stars 

^ . hosting a closely orbiting Jovian planet. The significance of this difference in 

metallicity arises from a homogeneous analysis of samples of FGK dwarfs which 
do not include the cooler and more problematic M dwarfs. This result would in- 
dicate that there is a possible link between planet mass and metallicity such that 
metallicity plays a role in setting the mass of the most massive planet. Further 
confirmation, however, must await larger samples. 
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Introduction 



More than 380 stars with planets have been discovered to date, half of which were 
detected in the past three years. Most of the extrasolar planets have been discovered via 
radial- velocity measurements of the reflex motions of the planet-hosting star and such surveys 
are biased to detect preferentially the largest and most closely orbiting planets. Within an 
ever increasing sample size, one statistically significant property has been confirmed for 
these objects: the average metallicity of the solar-like stars known to have giant planets (i.e. 
those planets close to the mass of Jupiter or larger) is h igher when compared to field F, G 
and K dwarfs no t known to host giant planets (see, e.g.. lGonzalezlll 997; 



Laws et al.ll2003l ; ISantos et al.ll2005l ; iFischer &: Valentil l2005t iBond et al. 



Santos et al. 



2006 



2001 



Sousa et al. 



20081 ). This difference is attributed to two possible scenarios: primordial enrichment and 
pollution. At present, the former seems to best account for the metal- rich nature of the 
planet hosting stars, since the probability of find i ng a planet is a steeply rising function 



of the stellar metallicity (e.g., ISantos et al.l l2004t Fischer fc Valentil 120051 ). However, the 



pollution hypothesis cannot be discarded, as contradictory conclusions have been found by 
several studies which attempted to unveil other chemical peculiarities in planet hosting stars 
(for a comprehensive review, see lGonzalezll2006l ; lUdry fc Santos! 120071 ) . 



In addition to the population of rather metal rich stars hosting Jovian-mass planets, 
there is a growing number of known systems with considerably lower-mass planetary com- 
panions. The range of planetary masses now includes objects with minimum masses of 
only about M p sini ~ 4M e , with many systems containing "Neptunian-mass" planets, with 
M p sin i < 25M e . It is of interest to investigate whether the trend for Jovian-mass planets 
to have a metal-rich stellar parent continues to wards systems wit h lower- mass planet s that 
do not conta i n the large Jovian-mass planets. lUdry et al.l ( 120061 ). ISousa et al.l (120081 ). and 
Mayor et al.l ( 120091 ) suggest that stars which have as their most massive planets Neptunian- 
mass objects may not be metal rich; however, the number of such systems which have been 
studied is just a few. The list of stars with Neptunian-mass planets continues to grow and 
these objects will help to probe the possibility of a stellar-metallicity planet-mass connection. 

The observed variety of exo-planetary masses and orbital separations, along with ev- 
idence of planetary migration and its possible influence on proto-planetary disk - stellar 
interactions suggests that it is of importance to determine chemical abundance distributions 
in different populations of exo-planet host stars. The search for subtle patterns in the abun- 
dances of stars with and without planets that may reveal details of planetary formation or 



1 Based on observations made with the 2.2 m telescope at the European Southern Observatory (La Silla, 
Chile), under the agreement ESO-Observatorio Nacional/MCT. 
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planetary system architecture is based ideally on a homogeneous and self-consistent analysis. 
If all samples are observed with the same instrumental setup and analyzed with a consis- 
tent methodology, systematic effects are more likely to be avoided. This study sets forth 
a homogeneous determination of stellar parameters and metallicities for a large sample of 
stars with planets, including a few stars hosting only Neptunian mass planets, as well as a 
control sample comprised of field stars not known to host giant planets. §2] describes the ob- 
servational data, sample selection criteria, and data reduction. The determination of stellar 
parameters, effective temperatures, surface gravities and metallicities, including the adopted 
iron line list, are presented in §[3j In §U results from this work are compared with those 
from the literature and discussed in light of various planet-metallicity correlations. Included 
in this discussion is an investigation into whether metallicity plays a role in determining the 
mass of the most massive planet in an exo-planetary system. Finally, concluding remarks are 
presented in §SJ The derivation of the elemental abundances other than Fe will be treated 
in subsequent papers. 



Observations and Data Reduction 



2.1. Sample Selection 



Planet Hosting Stars: The sample of main-sequence stars with planets analyzed in 
this study contains 117 targets. The target list was compiled using the Extrasolar Planet 
Encyclopaedia^ and updated with newly discovered systems until August 2008. We selected 
all planet hosting stars having 5 < +26° and V < 12. The declination limit was imposed 
by object observability at La Silla Observatory in Chile, and the limiting magnitude was 
set in order to keep exposure times needed to achieve the desired S/N relatively short. 
Sever al stars in this sample have been previously analyzed in recent studies of planet hosting 



stars (ILaws et al.l 



Luck fe Heiteril200d : 



200 



3; 



Santos et al. 



Bond et al.l 



2004 



200 



6: 



2005 



' 



Sousa et al 



Takeda et al. 2005; Valenti fe Fischer 2005 



2008 



;). We note that 16 planet hosting 



stars in this sample are not included in these previous abundance studies. 

Control Sample: A control sample of main-sequence disk stars whic h are not known to 
host g iant planets was selected from the list of nearby F, G and K stars in lFischer fe Valenti 
(120051 ) which has been targeted in the planet search programs conducted at the Keck Ob- 
servatory, Lick Observatory and the Anglo- Australian Telescope. That study identified 850 
stars for which there are enough observations to securely detect the presence of companions 



2 Available at http://exoplanet.eu 
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with velocity amplitudes K > 30 m s" 1 and orbital periods shorter than 4 yr. From the 
subsample of stars with non-detections of giant planets, we eliminated stars with [M/H] < 
-1.0; v sin i > 10 km s~ x (typical rotational velocities are much lower for solar-type stars) 
and 5 > +26°. In addition, any stars which were found subsequently to host giant planets 
(the only case being HD 16417) were obviously removed from the list. Binaries, as well as 
targets having one sing le spectrum analyzed were also excluded (according to Table 8 in 



Valenti fc Fischer 



Li (JGhezzi et al. 



2005). HD 36435 was added to the list as it was previously analyzed for 
20091 ). The final sample of comparison stars in this study has 145 tar- 
gets. A list with all targets analyzed, planet hosting stars as well as control sample stars, is 
presented in Table [TJ 



2.2. Observations and Data Reduction 

High-resolutio n spectra were obta ined with the Fiber-fed Extended Range Optical Spec- 



trograph fFEROS: lKaufer et al.lll999f ) attached to the MPG/ESO-2.20m telescope (La Silla, 
Chile). The detector was a 2k X 4k EEV CCD with 15 /xm pixels. This instrumental setup 
produces spectra with almost complete spectral coverage from 3,560 to 9,200 A (over 39 
echelle orders) and at a nominal resolution R = A/ A A ~ 48,000. The observations were 
conducted during 6 observing runs between April 2007 and August 200qj- A solar spectrum 
of the afternoon sky (T exp = 2x1 20s) was taken before each observing night. A detailed log 
of the observations, including V magnitudes, observation dates, total integration times and 
the resulting S/N per resolution element, is found in Table [TJ 

The spectra were reduced with the FEROS Data Reduction System (DRSJj The data 
reduction followed standard procedures. An average flat-field image was used in order to 
define the positions of the echelle orders. The background (bias level and scattered light) 
was subtracted from the images. The bias level was determined from the overscan region 
of the CCD and the scattered light was measured in the interorder space and in the region 
between the two fibers. The extraction of the echelle orders was done with a standard 
algorithm that also finds and removes cosmic rays. All extracted images were divided by 
the average flat-field in order to remove pixel-to-pixel variations and they were corrected 
for the blaze function. The flat-fielded spectra were wavelength calibrated using ThArNe 
and/or ThAr+Ne calibration frames. The calibrated spectra were rebinned in constant steps 
of wavelength and a barycentric correction was applied. Finally, the reduced spectra were 



3 Under the agreement ESO-Observatorio Nacional/MCT 

4 Available at http://www.eso.org/sci/facilities/lasilla/instruments/feros/tools/DRS/index.litml 
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corrected for radial velocity shifts by comparing the observed wavelengths of some isolated 
and moderately str ong iron lines with their rest wavelenghts taken from Vienna Atomic Line 
Databasfl (VALD; JKupka et al.lll999h . 



3. Analysis 



3.1. The Fe Line List 



The line lists for Fe I and F e II were compiled from the line sample in ISousa et al.l (120081 ) 
and iMelendez fe Barbuyl (120091) . The initial line list contained over 100 iron lines but using 
both the Solar Flux Atlas (IKurucz et al.l 119841 ) as well as the solar spectrum taken with 
FEROS on 20 August 2008, and from results of test calculations with a variety of (//-values 
from the literature, we selected 27 Fe I and 12 Fe II suitable lines which were unblended and 
of intermediate strength (equivalent widths less than 90 mA, in order to limit the effects of 
damping on the abundance determination). The final line list adopted in this Fe abundance 
analysis is presented in Table [2j The wavelengths and lower excitation potentials (LEP) of 
the Fe transitions were ta ken from VALD. T h e qf- values for Fe I tr a nsitio ns we re taken from: 



Blackwell et al. 



in 



(1982a. 



bl. ll984J. J1986J Il995h . bard. Kock fe Kockl Jl99lh and bard fe Kock" 



( 19941). and O'Brian et al.l ( 19911 ) . The gf -values from these stud ies were carefully comp ared 



Lambert et al.l ( 119961 ) and found to be in excellent agreement. iFulbright et al.l (120061 ) also 



argue that differences between the gf- values of these three groups are comparable to random 
uncertainties in the measurements. Corrections to the log gf scales from these different 
sources were therefore deemed not necessary in this study; whenever a transition had more 
than one gf- value available, an average value was ad opted. The qf- values for the Fe II lines 
in this study were taken from the critical analysis of lLambert et al.l (119961 ). 



3.2. Equivalent Width Measurements 



The code ARESn ( ISousa et al.l 120071 ) was used in order to measure equivalent widths 
(EWs) of sample Fe lines automatically. Briefly, this program first fits a polynomial to the 
local continuum in a spectral region defined by the user. It then determines which lines 
inside the given interval can be fit by a Gaussian profile. Finally, it computes the equivalent 
width(s) for the line(s) of interest assuming a Gaussian profile. More details about the 



5 Available at |http://ams.astro.univie.ac.at/~ vald/ 
6 http://www. astro. up.pt/ sousasag/ares/ 
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ARES code can be found in ISousa et al.l (120071 ) and ISousa et al.l ( 20081 ) . The equivalent 



widths measured for all program stars can be found in iGhezzil (J2010l ). 



Possible systematic effects in the automatic ARES equivalent width measurements were 
investigated here. We measured (using the splot task of IRAF) equivalent widths of 75 
Fe I lines and 22 Fe II lines in 6 stars which were selected to bracket the range in effective 
temperature, metallicity and spectrum quality (signal-to-noise ratio) of our sample as well 
as the Sun. A comparison between manual and automatic equivalent width measurements 
for 638 lines is shown in Figure [TJ The mean difference between the two sets of equivalent 
widths is (EW ares— ^^ Manual) = — 0.53±2.10 mA. Also, the following linear fit is obtained: 
EWiBEs = (1-005 ± 0.002)EW M arma« + (-0.77 ± 0.16). The correlation coefficient and the 
standard deviation are R = 0.99627 and a = 2.09 mA, respectively. 

The exercise above indicates that the equivalent widths which were measured automat- 
ically using the ARES code are consistent with our measurements, although there is a slight 
trend of ARES equivalent widths being marginally smaller th an the ones measur ed manually 
in this study. This result is in line with what was found in ISousa et al.l (120071 ) . Although 
we find an overall good agreement between manual and automatic equivalent widths, it is 
important to carefully check the results because ARES does not make quality assessments 
of the measurements it outputs. For instance, in the tests decribed above we note that 
there were 10 lines with obviously erroneous equivalent width measurements, which were 
discarded. In this study we estimate an uncertainty of ~2 mA as the typical uncertainty in 
the equivalent width measurements. Differences in equivalent widths of ± 2 mA are about 
what is expected given the resolution, sampling, and S/N of the spectra and no significant 
systematic effects are found between ARES and manual measurements. 



3.3. Derivation of Stellar Parameters and Iron Abundances 



Stellar parameters (T e ff, logg and £) and metallicities ([Fe/H]) were derived homo- 
geneously and following standard spectroscopic methods which are based on requirements 
of excitation and ionization equilibria. This abundance analysis wa s done in Loc al Ther- 
modynamic Equilibrium (LTE) using the 2002 version of MOOGLj ( Snedenl Il973l ). In all 
calculations van der Waals constants were multiplied by an enhancement factor of 2.0 
( JHolweger et al.lll99ll ). The model atmospheres in this study were interpolated from the 



7 Available at http://verdi.as.utexas.edu/moog.html 
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Fig. 1. — A comparison between manual (using IRAF splot) and automatic measurements 
of equivalent widths (using ARES code) for a sample of Fe I and Fe II lines in six target 
stars and the Sun. The solid line represents perfect agreement. 



ODFNEW grid of ATLAS9 model^l flCastelli fe Kurucz!l2004f ) 



Effective temperatures and microturbulent velocities were iterated until the slopes of 
A(Fe Ijj versus excitation potential, Xi an d A(Fe I) versus reduced equivalent width, 
log(EW/X), were respectively zero (excitation equilibrium). Only lines with \og(EW/X) < 
—5.00 (this limit was changed to larger values for cooler stars) were used in the first iteration, 
in order to decouple the T e ff an d £ determinations. Surface gravities were iterated until Fe I 
and Fe II returned the same mean abundances (ionization equilibrium). At the end of the 
iterative process, a consistent set of values of T e ff, log g and microturbulent velocity as well 
as the mean Fe I (= Fe II) abundance is obtained for the star. This procedure was adopted 
for all stars in our sample except for seven targets which had lower metallicities and solar 
temperatures (namely HD 6434, HD 51929, HD 80913, HD 114762, HD 153075, HD 155918 
and HD 199288). In these cases, the microturbulent velocities were kept at a fixed value 
because there were no lines strong enough in order to anchor the iteration of this parameter. 
As an example, Figure [2] shows the final iterated plots of A(Fe I) versus x (top panel) and 
A(Fe I) versus \og(EW/\) (bottom panel) for target HD 2039. 

An Automated Analysis: Due to the large number of stars in our sample and Fe lines 
included in the analysis, BASH and FORTRAN codes were built in order to automate the 
whole iterative process described above. In summary, the code starts with automatic equiv- 
alent width measurements using ARES (Section 13.21) and iterates to a final set of consistent 
values of effective temperature, surface gravity, microturbulence and Fe abundance (both 
from Fe I and Fe II). With the development of an automatic procedure it is now possible to 
analyze the entire sample of over 300 stars studied here in a few days without interventions. 

In order to further test our line list (Table [2]) and analysis method, the solar spectrum 
(observed with FEROS spectrograph on August 20, 2008) was analyzed in a similar manner, 
with automatic measurements of equivalent widths (the measured solar equivalent widths are 
found in the last column of Table [2]). Solar abundances A(Fe I) = 7.43 ± 0.07 and A(Fe II) 
= 7.44 ± 0.05 as well as a microturbulence £ = 1.00 km s _1 were derived using a Kurucz 
ODFNEW model atmosphere with T e g = 5777 K, log g = 4.44, a model turbulence of £ = 2.0 
k m s -1 , and 1/Hr> = 1 .25. T his solar Fe abu n dance is in excellent agreement with the results 
of lReddvetaD J2003I ) and iFulbright et all ( J2006h (A(Fe)=7.45), which use ^/-values from 



the same sources as here. The derived solar iron abundance also compares well (within the 
unce rtainties) with recent solar abundance determ inations for 3D hyd rodynamical models 
from lAsplund et all J2OO9I . A(Fe)= 7.50 ± 0.04) and lCaffau et~aD feoiol . A(Fe)=7.52 ± 0.06). 



8 Available at http://kurucz.harvard.edu/ 
9 A(Fe I) = log [N(Fe I)/N(H)] + 12 
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Fig. 2. — The spectroscopic determination of effective temperature and microturbulent ve- 
locity for HD 2039 obtained from zero slopes in the runs of Fe I abundances with excitation 
potential of the transitions (top panel) and reduced equivalent widths (bottom panel). The 
solid line represents the mean iron abundance and the dashed lines represent the 1-cr of the 
distribution. The top panel shows only those lines with log(EW/X) < —5.00 which were 
used in the T e ff iteration. Fe I and Fe II abundances are consistent and the slopes are zero 
for A(Fe) = 7.73. 
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Final values of effective temperatures, surface gravities and microturbuleny velocities 
for all stars are presented in Table |3j The metallicities [Fe/H] (listed in the last column 
in Table E]) were calculated for a solar abundance A(Fe) = 7.43 (as derived here). The 
a- values listed in columns 6 and 8 correspond to the standard deviations of the final mean 
abundances of Fe I and Fe II. The number of Fe I and Fe II lines considered for each star are 
listed in columns 7 and 9, respectively. 

As a comparison, photometric effective tem peratures were also deriv ed using the T e ff 
versus V — K calibration recently published by ICasagrande et al.l (120101 ). The V and K s 
magnitudes for the ta rget stars were taken, respectively, from The Hi pparcos and Tych o 
Catalogues JESAJIl997h and the 2MASS All-Sky Catalog of Point Sources JCutri et al.ll2003h . 
A comparison between the spectroscopic and photometric temperatures shows that the latter 
are systematically higher (Figure[3]). The mean difference for all studied stars is AT e ff(spec— 
phot) = —63 ± 113 K, indicating reasonable agreement. 



Uncertainties in the parameters T e g- , log g, £ and [Fe/H] were estimated as in lGonzalez fe Vanture 



(119981 ) and can be seen in Table [3J We note that these are internal errors and that the real 
uncertainties might be somewhat larger. Departures from LTE were not considered in this 
study and these can affect the derived LTE abundances. Non-LTE effects are expected to 
be smaller for Fe II lines as Fe II (and not Fe I) is the dominant ionization stage in solar 
type stars. For Fe I, departures from LTE are larger and may be at the level of ~0.1 dex 
(JGehren et alihoOlal lbl). 



3.4. Evolutionary Parameters 

Stellar luminosities, masses, radii and ages define the evolutionary stages of stars and 
for this sample these were calculated in the following way. Absolute magnitudes My were 
determined using the classical formula: 



M 



v 



\/ + 5 + 51og7T-v4 



V- 



(1) 



As already mentioned, the app arent V m agnitudes (Table [T]) were taken from The Hip- 
parcos and Tycho Catalogues (JESAl 119971 ). As the uncertainties in this magnitude are 
not listed, we adopted th e errors in Vr, given the similarities between the two passbands 
19971). The p a rallax es tt and their uncertainties (Table HJ columns 2 and 



(Ivan Leeuwen et al. 



3) were taken from Ivan Leeuwen! (120071 ). Three stars (namely HD 70573, BD-10 3166 and 
WASP 2) were not present in the sources above, thus their V magnitudes and parallaxes come 
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Fig. 3. — A comparison between photometric temperatures derived with the V — K cali- 
bration in ICasagrande et all (120101 ) and spectroscopic temperatures derived in this study for 
planet hosting stars (red open circles) and control sample (blue open squares). The solid line 
represents the bisector. The two effective temperature scales are in reasonable agreement, 
however, there is a tendency for the photometric T e fjs to be higher than the spectroscopic 
ones at the high T e ff end and lower at the low T e fj end. 
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from the references in SIMBAD 10 ! and The Extrasolar Planet Encyclopaedia. The interstellar 
monochromatic extinction at V magnitude, Ay ( Table HI column 4), as well as its error were 
calculated usi n g the tables from lArenou et al.l (11992) and the code EXTINCT. FOR from 
Hakkila et al.l (J1997I ). Note that the lArenou et al.l ( 119921 ) model is accurate to distances 
within 1 kpc of the Sun. 

Absolute magnitudes were converted to bolometric magnitudes M boi by adding bolomet- 
ric corrections in the V band, BCy (Table HJ column 5), linearly interpolated from the grids 
of iGirardi et al.l (120021 ) for the atmospheric parameters given in Table [3J The luminosities 
were then calculated using the well known relation: 



log^ = -0.4(M feoZ 



Mboi,e), 



(2) 



where M bo i,Q = 4.77 (IGirardi et al.ll2002l ). The uncertainties in My, M bo i and log(L/L©) were 
derived considering that a(BCy) and a(M bo i tQ ) are zero. The luminosities and uncertainties 
are listed in Table H] (columns 6 and 7). 

Effective te mperatures and l umino sities were used to place the stars on a grid of Y 2 
isochrones from iDemarque et al.l (120041), t hus allowing for an age determination. An inter- 
polation code provided by the authord_| was used in order to obtain a set of isochrones 
ranging from 1.0 to 13.0 Gyr in age and from —0.70 dex to +0.40 dex in metallicity (with 
steps of 1.0 Gyr and 0.1 dex, respectively). As an example, we show the grid of isochrones for 
[Fe/H] = +0.30 dex in Figure HI (top panel). The locations of all targets stars are indicated. 
The age of the closest isochrone was attributed for each star. Given the uncertainties in 
T e ff, l°g(L/L ) and [Fe/H] and the proximity of the isochrones, an age interval was also 
estimated for each star together with a single age (see Tabled]). Also, some stars were located 
outside the grid and their ages were indicated as lower or upper limits. 

Stellar radii and spectroscopic masses M spec (as well as their uncerta inties) were derived 
using standard relations (see e.g. eqs. 5-8 from lValenti fc Fischerl 120051 ) . The mass e s wer e 
also calculated by placing the stars on a grid of evolutionary tracks from lYi et al.l (120031 ) . 
The mass of the closest track was attributed for each star. An interpolation code (provided 
by the authors 12 !) was used in order to obtain a set of tracks ranging from 0.5 to 2.0 M Q in 
mass and from —0.70 dex to +0.40 dex in metallicity (with steps of 0.1 M & and 0.1 dex, 
respectively). A typical uncertainty in mass of 0.1 M & was estimated by considering the 



10 http://simbad. u-strasbg.fr/simbad/ 

11 Available at http://www.astro.yale.edu/demarque/yyiso.html 



12 Available at http://www.astro.yale.edu/demarque/yystar.html 
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errors in T e ff, log(L/L Q ) and [Fe/H]. In some cases, however, this error was larger because 
of the location of the star on the grid or due to a larger uncertainty in the luminosity. Also, 
a few stars were located below the Zero Age Main Sequence (ZAMS) and their masses had 
to be estimated through an extrapolation. As an example, the grid of evolutionary tracks 
for [Fe/H] = +0.30 dex is shown in Figure H] (bottom panel). The stellar radii and masses 
(spectroscopic, M spec , and those derived with the evolutionary tracks, M trac k) can be found 
in Table gj 

As the spectroscopic masses have greater errors, we adopt the masses obtained with 
the grid of evolutionary tracks. Using those masses, the Hipparcos surface gravities can be 
calculated with the relation: 

log # = log # + log— --log-— + 41og — — , (3) 

Mq Lq 1 C ff j0 

where T e s,& = 5777 K and logg = 4.44. The uncertainty in these gravities was calculated 
considering a(T e ff jQ ) = and o"(logg©) = 0. The results are presented in Table |4] (columns 
14 and 15). In Figure [5] we show a comparison between the derived spectroscopic log g's 
(listed in Table [3]) with Hipparcos log g's (listed in Table HJ column 14). The line indicating 
perfect agreement is also shown as a solid line in the figure. The agreement between the two 
sets of log g's is good although we note that the Hipparcos gravities are typically found to be 
higher (by 0.06 dex on the average) than the spectroscopic values with a standard deviation 
of ±0.15 dex, which is of the order of the estimated uncertainties in the derived log g from 
the iron line analysis. 

In addition, masses, radii, ages and trigonometric gravities were also derived with Leo 
Girardi's web code P ARAM 13 !. which is based on a Bayesian parameter estimation method 



( Ida Silva et al.l 120061 ) . The mean differences between the results discussed above (This work 
- Girardi's code) are small and indicate good agreement: AM = 0.03 ± 0.05M Q (N=262), 
AR = 0.01 ± 0.06^0 (N=223), At = 0.37 ± 1.46 Gyr (N=211) and Alogg = 0.03 ± 0.05 
(N=262). 



4. Discussion 

4.1. Comparisons with Other Studies 

Several recent studies in the literature have derived stellar parameters and metallicities 
for samples of planet hosting stars. In the following we briefly summarize some of these works 
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Available at http://stev.oapd.inaf.it/cgi-bin/param 
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Fig. 4. — The location of sample planet hosting stars (red circles) and control disk sample 
stars (blue squares) in an HR diagram. The values of effective temperature and luminosity 
for the targets are from Tables [3] and EJ respectively. The top panel shows isochrones for 
ages varying between 1 and 13 Gyr and the bottom panel shows evolutionary tracks for 
mass tracks between 0.8 and 2 Mm . The grids of isochrones and evolut ionary tracks were 
calculated for [Fe/H] = +0.30 dex JPemarque et~al~1l2004! : Iyj et al.ll2003h 
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Fig. 5. — A comparison between Hipparcos and spectroscopic gravities for planet hosting 
stars and control sample. The agreement between the two sets of surface gravities is found 
to be good and systematic differences between the two independent scales are less than ~ 0.1 
dex on average. The solid line is the bisector and represents the perfect agreement between 
the two determinations. 
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and then compare their results of effective temperatures, surface gravities and metallicities 
with the ones obtained in this study. 



Laws et a.1.1 (120031) determined spectroscopic parameters for 30 stars with giant plan- 



ets and/or brown dwarf companions. Their analysis method is similar to this study, the 
differ ence being the li ne list and qf -values which were obtained from an inverted solar anal- 



ysis. iSantos et a.1.1 ( 12004 . [2005) did a spectroscopic analysis of a large sample of stars 
with and without p lanets (119 and 94 targets, respectively). Their method is very similar 
to the one used by lLaws et al.l (20031). with th e difference being the list of iron lines (the 
gf- values are also solar). iTakeda et al.l (120051 ) obtained stellar parameters for a set of 160 
mid-F thr ough early-K dwarfs/subgiant s. The difference with previous studies is the selected 
iron lines. IValenti &c Fischerl ((2005) derived stellar properties for 1040 nearby F, G and 
K stars observed as part of the Keck Observatory, Lick Observatory and Anglo- Australian 
Telescope (AAT) planet search programs. Their method was different; stellar parameters 
and abundances were determined from a direct comparison of observed and synthetic spectra 
across certain spectral intervals using the spectral modelling program, SM E. In addition, a 
fixed y alue of 0.85 km s _1 for the microturbulence was adopted for all stars. iLuck fc Heiter 
(2006J) derived atmospheric parameters for a sample of 216 nearby dwarf stars. They used 
the standard spectro scopic method, but with differential abundances relative to the Sun. 
Bond et al.l ( 120061 ) determined atmospheric parameters from 136 G-type stars from the 
AAT planet search program. In that study, photometric temperatures are obtained from 
B — V colours listed in the Hipparcos catalog, with discrete values of microturbulence set 
at 1.00, 1.25 and 1.50 km s _1 . The metallicities and gravities are d etermined by i t erating 
this parameters until the Fe I and Fe II abundances are the same. ISousa et al.l (J2008) 
derived spectroscopic parameters for all 451 solar- type stars from the HARPS Guaranteed 
Time Observations (GTO) "hig h-precision" sample. Their method closely resembles that 
from ISantos et al.l ( 120041 . 120051 ). except for a larger line list and the usage of automatic 
measurements of equivalent widths. 

A direct comparison of the stellar parameters and metallicities obtained for some stud- 
ied stars in our sample with results from other studies discussed above is possible given that 
there are several targets in common. Table [5] shows the average differences (in the sense 
This study - Literature study) computed for the effective temperatures ((AT e ff)), surface 
gravities ((Alogg)) and met allicities ((A \ Fe/H ]) ) obtained fo r all target s ta rs we have in 



(2005); Valenti fe Fischerl (l2005l\ : ILuck fe Heiterl (j2006h 



comm o n with the studies of lLaws et al.l (120031); ISantos et al.l (12004 120051); ITakeda et al. 



Bond et al.l ( 120061 ); ISousa et al. 



( 120081 ); ISantos et al.l ( 120041 . 120051 ). The number of stars compared in each case is found 



in Table [5] (column 5). Results from this simple and direct comparison are briefly summa- 
rized below. 
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Effective Temperatures: In general, there is not a significant offset in the effective 
temperature scale in this study in comparison with the other studies in Table [5] In particular, 
for 5 studies we find (AT f .ff) to be less than 15 K, which is a quite small systematic offset; 
the iLuck &: Heiterl (120061 ) study has a diff erence which is onl y slightly larger (~35 K). A 
comparison with the T e ff from results in iBond et all (120061 ) indicates, however, a more 
significant systematic difference of ~75 K. The standard deviations around the mean values 
are all ~100 K or less, which is in general agreement with the estimated uncertainties for 
the derived effective temperatures in this study. 

Surface Gravities: A direct comparison between the average s urface gravity valu e 
derived for selected targ e ts in this study with a v erage results from lLaws et all (|20 03) ; 
Santos et all (booi J2005h : IValenti fc Fischer! (b()05h : ILuck fc Heiterl (|2006h and ISousa et al. 
( 120081 ) indicates that there is a small offset (~ 0.08 dex - 0.12 dex) in the lo gg scales. An 
agreem ent (at the level of 0.05 dex or better) is found between our results and iTakeda et al. 
( 120051 ) and IBond et all ( 120061 ). The stardard deviations of the distributions around the av- 
erage differences in log g are in all cases less than 0.2 dex; in agreement with the estimated 
uncertainties in the derived surface gravities. 

Metallicities: In terms of average metallicity values, the iron abundances derived in 
this study compare very well with results obtained in other studies in the literature for 
stars in common. There is a slight tendency, however, for the metallicities here to be just 
marginally lower (0.03 dex or less) than the other studies; but such differences a re probably 
statist ically insignificant. Note, however, that the iron abundance results in IBond et al. 
( 120061 ) are on average 0.09 dex lower than ours. 



4.2. Metallicity Trends with Effective Temperature and Stellar Mass 

Given our sample of 262 stars which have been subjected to a homogeneous analysis 
it is possible to search for differences in the properties of stars with planets compared to 
those stars not known to harbor giant planets with periods less than about 4 years. Two 
key defining properties of stars are their effective temperatures and masses, which on the 
main-sequence are related to each other, such that increasing T e g maps into increasing mass, 
at least over the relatively limited range of metallicities explored in this sample. In order to 
isolate possible differences between the two samples that might be related to T cff or mass, 
stars having surface gravities with log g < 4.2 were excluded from comparisons in this section. 
The resultant sample consists of 79 stars with planets and 109 stars without planets. Figure 
|6] compares the properties T e g and derived evolutionary track mass for the two samples of 
stars, with mass plotted versus T e g. The main-sequence nature of these stars is obvious from 
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the figure, with no significant differences in the distribution along the main sequence of stars 
with and without planets. 



4-2.1. Effective Temperatures, Iron Abundances and Solid-Body Accretion 

The defined set of target stars with log g > 4.2 (or those very near to the main sequence) 
are now illustrated with their values of [Fe/H] plotted versus T e g in Figure [7J The top panel 
contains stars without planets and the bottom panel shows stars with planets. Since these 
are main-sequence stars, the effective temperature follows the stellar mass. No strong trends 
between [Fe/H] and T e g are apparent in either sample, with the lack of an increase of [Fe/H] 
with increasing effective temperature placing limits on the amount of solid-body accretion 
(material depleted in H and He, for example) that might have occurred in these stars, 
since the convective-zone mass of a main-sequence star is a strongly decreasing function of 
increasing T e g. 



This same accretion test was conducted by iPinsonneault et al.1 (1200 if ) on an early small 
sample (~30) of stars with planets and they found no trend of increasing [Fe/H] with T eff . 
Accretion of solid material could create a positive trend due to the significantly decreasing 
convective zone mass in main-sequence stars with increasing effective tempratures: e.g., the 
co nvective-zone mass decrea ses by about a factor of 50 in going from T c g = 5000 K to 6400 



K (IPinsonneault et al.1 120011 ) . Accretion of only a few Earth masses of solid material would 



increase the surfac e value of [Fe/H] by ~+0. 3 dex in a solar- metallicity star with T c g = 6400 



K (see Figure 2 in IPinsonneault et al.ll200lf ). No such trend is seen in Figure [7J suggesting 
that accretion of more than a few Earth masses of solid material is either rare, or such 
accreted material sinks rapidly out of the outer convection zone. 



4-2.2. Stellar Mass and Metallicity 

In addition to the comparison carried out above between [Fe/H] and T e e, it is also 
instructive to do a similar comparison with stellar mass (in this case using the evolutionary 
track masses); this comparison is shown in Figure [HI again, where stars without planets are 
plotted in the upper panel and stars with planets in the lower panel. The samples have been 
binned in mass intervals of 0.25 M , as represented by the error bars in the abscissa. The 
values of [Fe/H] plotted represent the mean value within that mass interval, with the error 
bars showing the standard deviations of [Fe/H] at that mass. In both samples, the values 
of [Fe/H] increase with increasing stellar mass. Such an increase was noted in the review by 
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Fig. 6. — The relation between evolutionary track mass and effective temperatures for stars 
with (red open circles) and without (blue open squares) planets. Stars in both samples 
populate a well-defined "main- sequence" , with significant overlap in mass and T e ff between 
the stars with giant planets and those without. 
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Fig. 7. — The trend between metallicities and effective temperatures for control sample stars 
(upper panel) and planet hosting stars (lower panel). There is a spread of metallicities at 
each T e ff, with no strong trends. Significant solid-body accretion would produce an upturn 
in the upper envelope of [Fe/H] with increasing T e ff'- such an effect is not observed in either 
sample. 
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Gonzalezl (120061 . his Figure 1) using the abundance results from iFischer fc Valentil (120051 ) 



A signature of solid-body accretion polluting the stellar convective envelopes would be 
an upturn in [Fe/H] with increasing stellar mass, since the convective envelope mass is a 
rapidly decreasing function of increasi ng stellar mass. At first glanc e, the increase in [Fe/H] 



with mass found here (and noted by iGonzalea 120061 based on the IFischer fc Valentil 12005 



results) might suggest that solid-body accretion has taken place. Two effects, however, 
indicate that this has not affected the overall metallicities. First, the slopes of A [Fe/H] /AM 
are identical in both the stars with planets and stars without planets. This slope is very 
ro ughly +0. 7 dex/ Mm and is similar to the slope that would be deduced from Figure 1 



m 



Gonzalezl ( 120061 ). The fact that all of the various samples of stars, with and without 



giant planets, exhibit similar behavior in metallicity with mass argues that pollution has not 
selectively altered the values of [Fe/H] in a significant way for the stars with giant planets. 
The second point to note is that a positive slope of [Fe/H] with stellar mass would result from 
an age-metallicity relation. Since more massive stars have shorter main-sequence lifetimes, 
they would be biased towards higher values of [Fe/H], while lower-mass stars would be a 
mixture of old and young stars, which would shift the overall distribution to lower average 
values of [Fe/H]. 

In summary, comparisons among iron abundances with effective temperatures and stel- 
lar masses in both samples of stars (with and without large planets, respectively) reveal that 
accretion of solid-body material does not affect significantly (~0. 1-0.2 dex) the overall bulk 
metallicity in either sample. This does not rule out smaller amounts of accretion, which 
might affect abundance ratios bet ween certain type s of elements (such as vo latile versus 



refractory species, as suggested by lSmith et al.ll200ll ; see iMelendez et al.ll2009l for an alter 



native interpretation). This question will be addressed in a later paper using the spectra 
from this dataset and analyzing a broad range of elements. 



4.3. Metallicity Distributions of Sample Stars 



As discussed in Section l4~Tl the metallicities ([Fe/H]) derived here are generally consistent 
(within the expected errors) with metallicities found in other studies of planet hosting stars in 
the literature. As the present study relies on a homogeneous and self-consistent analysis of a 
sample of 262 stars, having comparable numbers of planet hosting and comparison disk stars, 
it is possible to quantify differences in the metallicity distributions in these two populations. 
Figure [9] shows the metallicity distributions for stars with planets (solid line histogram) and 
comparison stars (dashed line histogram). There is an offset in the peak metallicity of the 
two histograms in the figure. The peak of the distrbution for stars with planets is located 
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Fig. 8. — Average metallicities versus binned evolutionary track mass for control sample 
stars (upper panel) and planet hosting stars (lower panel). Masses are binned in 0.25 M 
intervals with each point showing the mean value and standard deviation of [Fe/H] within 
that mass interval. The slopes of A [Fe/H] /AM are the same within their uncertainties, 
in each sample and are likely the result of the age-metallicity relation. Selective accretion 
of solid-body material in stars with giant, closely orbiting planets would result in different 
slopes between the two samples and this is not observed. 
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in the bin centered at [Fe/H] = +0.20 dex and the average metallicity of this sample is 
([Fe/H]) = +0.11 dex. For the comparison stars, the peak is on the bin centered at [Fe/H] 
= +0.10 dex and the average metallicity in this case is lower: ([Fe/H]) = —0.04 dex. Thus, 
there is an offset of 0.15 dex between the mean metallicities of the two samples. 

When comparing properties, such as metallicity, between samples of planet hosting stars 
and those without giant pla nets, it is worth noting some selection biases inherent in these 



samples. As summarized by iGonzalezI (120061 ) . Doppler surveys avoid young, chromospher- 



ically active stars (which also typically are fast rotators) and contain only small numbers 
of metal poor stars (with [Fe/H] < -0.5 dex) because these objects are rare in the solar 
neighbourhood. In particular, our control sample of disk stars would also suffer from such 
biases since it was selected in order to search for the presence of planets. 

As this offset between the peak values of the two histograms is of the order of, or just 
slightly higher than, the expected uncertainties in the derived iron abundances themselves, 
it is useful to perform a robust statistical test in order to further investigate whether the 
metallicity offset is meaningful. In this sense, we conducted a Two-Sample Kolmogorov- 
Smirnov test and found a probability P = 6.17 x 10~ 6 that the two samples are drawn from 
the same parent population. This low probability confirms the results previously found in 
the literature that the population of stars hosting giant planets is more metal-rich than the 
population of stars not known to harbour such planets. 

If a volume limited sample with a radius of 18 pc is defined here for comparison, the mean 
metallicity for the control sample disk stars (N=46) is now -0.11 dex and the offset relative 



to pla net-hosting stars becomes 0.22 dex, similar to the one found by iFischer fc Valenti 



( 120051 ) based on a much larger sample. The average metallici ty found for the volume-lim ited 



sample in this stu dy is also consistent with the results from ISantos et al.l ( 12004 . 120051 ) and 



Sousa et al.l ( 120081 ). The former study uses a comparison sample of 94 stars within 20 pc of 
the Sun and finds ([Fe/H]) = —0.11 dex. The latter work extends the comparison sample 
to 385 stars and the enclosed radius to 56 pc and finds ([Fe/H]) = —0.12 dex. 

The results from Sections 14.2. II and I4.2.2l indicate that solid body accretion has probably 
not altered surface values of [Fe/H] at the level of the offset in metallicity; the difference 
in [Fe/H] between the two samples suggests that intrinsic metallicity influences giant planet 
formation and migration. 
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Fig. 9. — Metallicity distributions obtained for planet hosting stars (red solid line) and a 
control sample of disk stars not known to host giant planets (blue dashed line). The peaks of 
the metallicity distributions are offset by 0.15 dex indicating that the sample of stars which 
host planets is typically more metal rich. 
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4-3.1. Metallicities of stars hosting Neptunian-mass planets 

The conclusions from the previous section favor the premise that metallicity plays a role 
in influencing the formation of giant planets, i.e., planets with roughly the mass of Jupiter. 
Within this model it is worthwhile exploring whether stellar metallicity also plays a role in 
the mass distribution of planetary systems. Such a comparison begins to probe how the 
underlying planetary system architecture might depend on metallicity. In this section, the 
metallicities of stars harboring only lower mass planets, i.e. those having Neptunian masses, 
are compared to systems containing the larger Jovian mass planets. 



Sousa et al.l ( 120081 ) presented preliminary results in which they found possible metallicity 



differences between stars hosting Jovian mass planets compared to those hosting Neptunian- 
mass planets (M p sini < 25M ffi ). The differences between the two metallicity distributions 
defined respectively by stars with Jovian-mass planets as opposed to Neptunian-mass planets 
indicated that these groups are not likely to belong to the same populations of stars. This 
conclusion, however, was based on a comparison of 63 Jovian hosting stars with a sample 
of 11 Neptunian hosting stars (those which contained at least one Neptune mass planet). 
Five of these 11 stars were M dwarfs and their metallicities were taken from the literature, 



while three were FG dwarfs analyzed by ISousa et al.l (120081 . HARPS GTO). The difference 
between the two metallicitiy distributions was a small mean offset of about 0.1 dex, with the 
Neptune hosting stars having a slightly lower mean metallicity. 

Because M dwarfs are more difficult to analyze spectroscopically, due to considerable 
line blending and blanketing from molecules, it is necessary to consider abundance uncer- 
tainties and systematics when using res ults from the s e complex stellar spectra . Recent 



abundance analys es of M dwarfs include IBonfils et al.l (120051 ) ; I Woolf fc Wallersteinl (120061) ; 
Bean et all (l2006h and I johnson fc Ap psl (l2009h . The studies bv IBonfils et~aD J2OO5I ) and 
Johnson fc Appsj (120091 ) point, fo r instance, to potentially large uncertainties in derived M 
dwarf abundances. For example, Ijohnson fc Appd (120091 ) find an average off set of ~ 0.3 



dex in [Fe/H] when their abundances are compared to the same M dwarfs from IBonfils et al. 



( 120051 . see Table E]). Such discrepancies suggest that, until results for the cooler M dwarfs 
are on firmer ground, it is prudent to investigate the effects of both including M dwarf 
metallicities in such comparisons, as well as excluding them. 

The sample studied here contains 9 systems which host at least one Neptunian mass 
planet, none of which are M dwarfs. This is the largest sample of stars hosting Neptunian 
size planets analyzed homogeneously for metallicities to date and this subsample can be 
directly compared to the Jupiter-like planet hosting stellar sample. The strength of such 
a comparison rests upon the fact that all stars have been analyzed homogeneously and 
are within a similar range of stellar parameters, so that systematic errors are not likely to 
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produce spurious differences in the metallicity distributions. The weakness is that the sample 
of Neptune-mass hosts has only a small number of stars. 

Figure [TD] (top panel) shows histograms representing metallicity distributions of two 
samples: those stars hosting at least one Jupiter-mass planet (N=112; black solid line) and 
stars hosting only Neptune-mass planets (N=5; represented by the dashed red line). There 
is a hint that stars with only Neptunian planets tend to be more metal poor compared to 
Jovian-planet hosting stars. The average metallicity of the Neptunian hosts is -0.08 dex, 
while the Jupiter host metallicity distribution has an average of +0.12 dex. If a Two- 
Sample Kolmogorov-Smirnov test is performed, we find a probability of 8% that the stars in 



the t wo samples belong to the same metallicity population (which agrees with ISousa et al. 



20081 ). This is a tantalizing result that suggests that metallicity may play a role not just in 
the formation of giant planets, but may also influence the distribution of planetary masses 
within exo-solar systems. This important question needs to be answered more definitively, 
but this will require larger samples. 

In order to extend the sample of stars with Neptunian mass planets (shown in the top 
panel of Figure [TU]) in this discussion, a list of stars with at least one Neptunian planet was 
compiled from The Extrasolar Planet Encyclopaedia and is presented in Table |6] along with 
the metallicity results from the different studies in the literature. In a similar analysis as 
was done for the sample of stars studied here (discussed above), Two-Sample Kolmogorov- 
Smirnov tests were done now with the inclusion of the literature sample using different 
permutations. The results from these tests are discussed below. 

1) Using the literature values of [Fe/H] for only F, G, and K dwarfs (no M dwarfs) a dif- 
ference in the mean [Fe/H] of +0.11 dex (in the sense of Jovian-mass hosts minus Neptunian- 
mass hosts) is found, with a probability of P=17% that the two samples were drawn from 
the same [Fe/H] populations (with N( Jovian- mass) =112 and N(Neptunian-mass)=15). The 
histogram showing the comparison of these two distributions is presented in the bottom 
panel of [TOl 



2) Using all literature values, with M dwarf abundances from iJohnson fc Appsj (120091 ) 
included, we find A[Fe/H]=+0.10 dex and P=17% (N(Jovian-mass=112 and N(Neptunian- 
mass)=19). 



3) When M-dwarf abundances from iBonfils et al.l (120051 . 120071 ) are used instead 



A[Fe/H]=+0.14 dex and P=5% (N(Jovian-mass)=112 and N (Neptunian-mass) =18). 

In the above exercise the stars with planets were divided into systems with Jovian-mass 
planets and those with Neptunian-mass planets, respectively. The metallicity comparison 
can also be carried out by dividing the sample into stars with at least one Neptunian-mass 
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Fig. 10. — A comparison between the metallicity distributions for planet hosting stars. Top 
panel: The black solid line histogram represents the sample of planet hosting stars containing 
at least one Jupiter planet. The red dashed line represents the sample of stars hosting only 
planets with Neptune-like masses. The metallicities are from this study. Bottom panel: 
Results for metallicities in Neptunian-mass hosting stars from the literature are added (red 
dashed line histogram). Given the large uncretainties in the M dwarf metallicities these 
stars have not been added to the sample. There is a hint that there is offset between the 
two distributions represented both in the top and bottom panels. 



-28- 



planet, regardless of whether there is also a Jovian mass planet in the system, and those 
systems with Jovian mass planets but no Neptunian-mass planets. The average metallicity 
of the sample of stars which host at least one Neptunian planet (N=9) is +0.09 dex, close to 
the value derived for stars hosting at least one Jupiter-mass planet (+0.12 dex; N=112). A 
Two-Sample Kolmogorov-Smirnov test reveals that the probability that these two samples 
belong to the same parent population is 91%. This comparison strengthens the idea that 
metallicity played a role in setting the mass of the most massive planet in a system. 

All of these various comparisons taken together suggest that lower values of [Fe/H] 
tend to produce lower-masses of the most massive planet within a planetary system. Such 
conclusions, however, should be viewed with caution since the number of stars harboring 
Neptune-size planets analyzed to date is still rather small, b ut in line with what would be ex- 



pected from models of plan et formatio n via core accretion (jlda fc Linll2004j ; iMordasini et al. 
2009al B see also review by lBossllioiol ). 



5. Conclusions 

We have determined stellar parameters for 117 main-sequence stars with planets discov- 
ered via radial velocity surveys and 145 comparison stars which have been found to exhibit 
nearly constant radial velocities and are not likely to host large, closely orbiting planets. 
The stellar parameters were derived from a classical spectroscopic analysis, using accurate 
laboratory gf- values of Fe lines and automatically measured equivalent widths, after critical 
evaluations of their quality. The values of effective temperature, surface gravity (as logg), 
microturbulent velocities, and iron abundances are in general good agreement with most of 
the values presented in a number of literature studies, but problems with a few individual 
stars may remain. 

Correlations between [Fe/H] and T e g in members of either sample are not found, which 
places stringent limits on the possible accretion of solid material (of less than a few Earth 
masses) onto the surfaces of these stars. A trend of increasing [Fe/H] with increasing stellar 
mass is found in both samples of stars, with the slope of A [Fe/H] /AM being the same 
for stars with giant planets and control sample stars. The same value of A [Fe/H] /AM for 
both samples rules out solid-body accretion, which leaves the underlying disk age-metallicity 
relation as the likely cause of the positive correlation of [Fe/H] with mass. 

It should be noted that the samples analyzed here were not selected based on any 
rigourous criteria, other than being segregated based on the presence of giant planets in 
one sample and the probable absence of such planets in the other. The list of stars without 
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planets preferentially included metal-rich stars, while some stars with planets were discovered 
by surveys that were based on high-metallicity as a criterion. 

A comparison of iron abundances between the stars with planets with the results from 
the control sample of disk stars confirms the results obtained by previous studies showing 
that planet-hosting stars exhibit larger metallicities than stars not harboring planets; the 
difference found in this sample is that stars with planets are shifted by +0.15 dex in [Fe/H] 
when compared to the stars without planets. 

The sample of stars with planets discussed here contains 9 stars which host at least 
one Neptunian- mass planet; of these 9 systems, 4 also contain a Jovian-mass planet, while 
5 contain only Neptunian-sized planets and no Jovian-mass closely orbiting planets. A 
statistical test of the iron abundances indicates that there is probably a real difference 
between the metallicity distributions of stars which contain only smaller Nep tunian-sized 



plane ts in comparison with stars hosting the larger Jovian- mass planets (see also lSousa et al. 



2008|). 



Although it should be recognized that the sample sizes are still small, there seems 
to be early indications that metallicity plays an important role in setting the mass of the 
most massive planet. It is also important to note that such a conclusion obtained here 
is based on stars which have all been analyzed homogeneously and which also have similar 
stellar parameters, therefore avoiding the large uncertainties still hampering the analyses and 
derived metallicities for the cooler M dwarfs. The same statistical test applied to metallicity 
results in the literature obtained for all FGK stars which host only Neptunian-sized planets 
allows for a larger sample and overall corroborates the results obtained with our sample. 

We would like to dedicate this paper to the co-author Chico Araiijo who was our friend, 
adviser and colleague and passed away in 2009. We acknowledge the financial support of 
CNPq. Research presented here was supported in-part by NASA grant NNH08AJ581. Luan 
Ghezzi thanks Oliver Schiitz for his valuable help in using the FEROS DRS package, and 
Claudio Bastos for conducting the observing run in February 2008. We thank the anonymous 
referee for useful comments that helped improving the paper. 
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Table 1. Log of Observations 



Star 


V 


Observation 


T 

■*- exp 


S/N 






Date 


(s) 


(~ 6700 A) 


Planet Hosting Stars 






HD 142 


5.70 


2007 Aug 30 


200 


367 


HD 1237 


6.59 


2007 Aug 30 


480 


466 


HD 2039 


9.00 


2007 Aug 28 


3000 


362 


HD 2638 


9.44 


2007 Aug 29 


3000 


276 


HD 3651 


5.88 


2007 Aug 29 


200 


356 


Control Sample 








HD 1581 


4.23 


2008 Aug 20 


30 


216 


HD 1835 


6.39 


2007 Oct 02 


200 


388 


HD 3823 


5.89 


2007 Oct 02 


200 


452 


HD 4628 


5.74 


2008 Aug 20 


200 


375 


HD 7199 


8.06 


2008 Aug 19 


1200 


348 



Note. - - Table [T] is published in its entirety in the 
eletronic edition of the Astrophysical Journal. A por- 
tion is show here for guidance regarding its form and 
content. 
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Table 2. Selected Fe lines and Measured Equivalent Widths for the Sun. 



A Ion LEP log gf W A 

(A) (eV) (dex) (mA) 



4779.439 


Fel 


3.415 


-2.020 


40.4 


4788.751 


Fel 


3.237 


-1.763 


63.3 


4802.875 


Fel 


3.695 


-1.514 


58.3 


4962.572 


Fel 


4.178 


-1.182 


52.7 


5054.642 


Fel 


3.640 


-1.921 


39.6 


5247.049 


Fel 


0.087 


-4.946 


65.8 


5379.574 


Fel 


3.694 


-1.514 


59.8 


5618.631 


Fel 


4.209 


-1.276 


48.4 


5741.846 


Fel 


4.256 


-1.670 


30.8 


5775.081 


Fel 


4.220 


-1.298 


57.9 


5778.453 


Fel 


2.588 


-3.430 


21.3 


5855.076 


Fel 


4.608 


-1.478 


21.4 


5916.247 


Fel 


2.453 


-2.994 


55.0 


5956.692 


Fel 


0.859 


-4.605 


51.3 


6120.246 


Fel 


0.915 


-5.970 


5.2 


6151.617 


Fel 


2.176 


-3.299 


48.7 


6173.334 


Fel 


2.223 


-2.880 


67.2 


6200.313 


Fel 


2.608 


-2.437 


72.5 


6219.279 


Fel 


2.198 


-2.433 


89.8 


6240.645 


Fel 


2.223 


-3.170 


48.0 


6265.131 


Fel 


2.176 


-2.550 


84.2 


6380.743 


Fel 


4.186 


-1.376 


51.8 


6593.870 


Fel 


2.433 


-2.422 


84.5 


6699.141 


Fel 


4.593 


-2.101 


7.7 


6739.520 


Fel 


1.557 


-4.794 


11.7 


6750.150 


Fel 


2.424 


-2.621 


73.2 


6858.145 


Fel 


4.607 


-0.930 


51.1 


4993.358 


Fell 


2.807 


-3.670 


37.6 


5132.669 


Fell 


2.807 


-4.000 


26.9 


5284.109 


Fell 


2.891 


-3.010 


62.6 
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Table 2 — Continued 



A Ion LEP log gf W A 

(A) (eV) (dex) (mA) 



5325.553 


Fell 


3.221 


-3.170 


39.2 


5414.073 


Fell 


3.221 


-3.620 


27.6 


5425.257 


Fell 


3.199 


-3.210 


41.5 


5991.376 


Fell 


3.153 


-3.560 


30.4 


6084.111 


Fell 


3.199 


-3.800 


20.5 


6149.258 


Fell 


3.889 


-2.720 


35.0 


6369.462 


Fell 


2.891 


-4.190 


19.9 


6416.919 


Fell 


3.892 


-2.680 


38.5 


6432.680 


Fell 


2.891 


-3.580 


40.4 



Table 3. Atmospheric Parameters. 



Star 


T e// 


log g 


£ 


A(Fe) 


a 


N 


a 


N 


[Fe/H] 




(K) 




(km s _1 ) 




(Fel) 


(Fel) 


(Fe II) 


(Fe II) 




Planet Hosting Stars 


















HD 142 


6338 ± 46 


4.34 ± 0.14 


2.27 ± 0.08 


7.46 


0.09 


23 


0.07 


10 


0.03 ± 0.04 


HD 1237 


5572 ± 40 


4.58 ± 0.09 


1.34 ± 0.04 


7.55 


0.09 


27 


0.05 


9 


0.12 ± 0.04 


HD 2039 


5934 ± 36 


4.30 ± 0.13 


1.26 ± 0.04 


7.73 


0.08 


27 


0.05 


12 


0.30 ± 0.03 


HD 2638 


5236 ± 70 


4.38 ± 0.19 


0.86 ± 0.04 


7.65 


0.07 


26 


0.07 


9 


0.22 ± 0.05 


HD 3651 


5252 ± 65 


4.32 ± 0.18 


0.81 ± 0.02 


7.62 


0.08 


27 


0.05 


10 


0.19 ± 0.03 


Control Sample 


















HD 1581 


5908 ± 31 


4.26 ± 0.13 


1.17 ± 0.04 


7.23 


0.08 


25 


0.08 


12 


-0.20 ± 0.03 


HD 1835 


5829 ± 41 


4.39 ± 0.16 


1.24 ± 0.04 


7.65 


0.07 


22 


0.05 


10 


0.22 ± 0.03 


HD 3823 


6012 ± 31 


4.18 ± 0.08 


1.92 ± 0.10 


7.08 


0.07 


25 


0.05 


11 


-0.35 ± 0.02 


HD 4628 


5055 ± 40 


4.33 ± 0.19 


0.88 ± 0.04 


7.14 


0.08 


25 


0.06 


5 


-0.29 ± 0.02 


HD 7199 


5349 ± 65 


4.09 ± 0.19 


1.04 ± 0.04 


7.74 


0.09 


26 


0.05 


10 


0.31 ± 0.05 



Note. - - Table [3] is published in its entirety in the eletronic edition of the Astrophysical Journal. A portion 
is show here for guidance regarding its form and content. 
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Table 4. Evolutionary Parameters. 



Star 


7T 


o~ n 


A v 


BC V 


log(L/L Q ) 


CT log(L/L ) 


R 


O-R 


■Lvlspec 


a(M apec ) 


Mtrack 


°~(M tr ack) 


log gHip P 


o-(log gHipp) 


Age 


AAge 




(mas) 


(mas) 










(Re) 


(Re) 


(M ) 


(M Q ) 


(M Q ) 


(M ) 






(Gyr) 


(Gyr) 


Planet Hosting Stars 
































HD 142 


38.89 


0.37 


0.05 


0.017 


0.462 


0.061 


1.41 


0.11 


1.59 


0.77 


1.25 


0.10 


4.24 


0.08 


2.5 


2.0-3.5 


HD 1237 


57.15 


0.31 


0.04 


-0.075 


-0.196 


0.060 


0.86 


0.07 


1.02 


0.49 


1.00 


0.10 


4.57 


0.08 


< 1.0 


0.0-3.0 


HD 2039 


9.75 


0.95 


0.11 


-0.004 


0.376 


0.104 


1.46 


0.18 


1.55 


0.81 


1.25 


0.10 


4.21 


0.11 


3.5 


3.0-4.5 


HD 2638 


20.03 


1.49 


0.10 


-0.159 


-0.368 


0.089 


0.80 


0.09 


0.55 


0.28 


0.90 


0.10 


4.59 


0.11 


< 1.0 


0.0-4.0 


HD 3651 


90.42 


0.32 


0.02 


-0.153 


-0.287 


0.060 


0.87 


0.07 


0.57 


0.28 


0.90 


0.10 


4.51 


0.08 


5.0 


0.0-12.0 


Control Sample 
































HD 1581 


116.46 


0.16 


0.01 


-0.037 


0.102 


0.060 


1.08 


0.08 


0.76 


0.37 


1.00 


0.10 


4.38 


0.08 


6.5 


4.0-9.0 


HD 1835 


47.93 


0.53 


0.08 


-0.025 


0.033 


0.061 


1.02 


0.08 


0.93 


0.45 


1.10 


0.10 


4.47 


0.08 


1.0 


0.0-3.5 


HD 3823 


40.07 


0.34 


0.04 


-0.035 


0.376 


0.061 


1.42 


0.11 


1.11 


0.54 


1.00 


0.10 


4.13 


0.08 


7.5 


6.5-9.0 


HD 4628 


134.14 


0.51 


0.01 


-0.226 


-0.549 


0.060 


0.69 


0.06 


0.37 


0.18 


0.70 


0.10 


4.60 


0.09 


9.0 


0.0-14.0 


HD 7199 


28.33 


0.57 


0.10 


-0.122 


-0.132 


0.063 


1.00 


0.08 


0.45 


0.22 


0.95 


0.10 


4.41 


0.08 


9.0 


5.0-13.0 



CO 

oo 



Note. 



Table|4]is published in its entirety in the eletronic edition of the Astrophysical Journal. A portion is show here for guidance regarding its form and content. 
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Table 5. Comparison with other results in the literature. 



Study 



(AT e// )(K) (A logs) (A[Fe/H]> N Stors 



Laws et al. f2003) 
Santos et al. f2004 . 2005) 
Takeda et al. (20051 
Valenti fc Fischer f2005) 
Luck fc Heiter f 20061 
Bond et al. (20061 
SousaeLaL C2008 ) 



-5±74 


-0.10±0.15 


-0.03±0.05 


23 


-2±72 


-0.08±0.13 


-0.02±0.06 


113 


-8±65 


-0.05±0.16 


-0.03±0.07 


35 


10±65 


-0.11±0.14 


-0.01±0.06 


223 


-32±84 


-0.11±0.15 


-0.02±0.07 


56 


74±113 


0.01±0.19 


0.09±0.09 


90 


-14±61 


-0.11±0.11 


0.00±0.06 


119 



Note. - - A = This study - Literature Study 
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Table 6. Neptunian-mass Planet Hosts 



Star 


Mp sin i 

(M e ) 


Jupiter 


[Fe/H] 


Reference 

[Fe/H] 


Results from 


This Work 








HD 4308 


12.87 


no 


-0.31 




HD 16417 


21.93 


no 


0.14 




HD 40307 


4.20 


no 


-0.35 




HD 47186 


22.78 


yes 


0.21 




HD 69830 


10.49 


no 


0.00 




HD 125612 


21.29 


yes 


0.25 




HD 160691 


10.56 


yes 


0.23 




HD 181433 


7.56 


yes 


0.46 




HD 219828 


20.98 


no 


0.14 




Literature Results 








HD 7924 


9.22 


no 


-0.15 


Howard et al. (2009) 


HD 1461 


7.60 


no 


0.18 
0.21 
0.19 
0.19 


Valenti & Fischer (2005) 
Luck & Heiter (2006) 
Sousa et al. (2008) 
Average 


CoRoT-7 


4.80 


no 


0.05 


Leeer et al. (2009) 


55 Cnc 


7.63 


yes 


0.33 


Santos et al. (2004) 


BD-082823 


14.30 


yes 


-0.07 


Hebrard et al. (2010) 


HD 90156 


17.48 


no 


-0.24 


Encyclopaedia 


61 Vir 


5.09 


no 


0.01 
0.05 
0.11 
-0.02 
0.04 


Santos et al. (2004. 2005) 
Takeda et al. (2005) 
Valenti & Fischer (2005) 
Sousa et al. (2008) 
Average 


HD 125595 


14.30 


no 


0.02 


Encyclopaedia 


HD 156668 


4.16 


no 


-0.07 


Mishenina et al. (2008) 


Kepler-4 


24.47 


no 


0.17 


Borucki et al. (2010) 


HD 179079 


25.43 


no 


0.25 


Valenti et al. (2009) 
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Table 6— 


-Continued 


Star 


Mp sin i 

(M e ) 


Jupiter 


[Fe/H] 


Reference 

[Fe/H] 


HAT-P-11 


25.74 


no 


0.31 


Bakos et al. (20101 


HD 190360 


18.12 


yes 


0.24 


Sousa et al. (2008) 


HD 215497 


5.40 


yes 


0.23 


Lo Curto et al. (2010) 


Literature Results for . 


M Stars 






HD 285968 


8.42 


no 


-0.10 
0.18 
0.04 


Endl et al. (2008) 
Johnson & Adds (2009) 
Average 


GJ436 


22.88 


no 


0.02 
-0.32 

0.25 
-0.02 


Bonfils et al. (2005) 
Bean et al. (2006) 
Johnson & Adds (2009) 
Average 


G1581 


1.94 


no 


-0.25 
-0.33 
-0.10 
-0.23 


Bonfils et al. (2005) 
Bean et al. (2006) 
Johnson & Adds (2009) 
Average 


GJ674 


11.76 


no 


-0.28 
-0.11 
-0.20 


Bonfils et al. (2007) 
Johnson & Adds (2009) 
Average 


Gliese 876 


6.36 


yes 


-0.03 

-0.12 

0.37 

0.07 


Bonfils et al. (2005) 
Bean et al. (2006) 
Johnson & Adds (2009) 
Average 



